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Abstract  Progress  is  reported  in  development,  implemen¬ 
tation,  and  application  of  a  spectral  method  for  ab  initio 
studies  of  the  electronic  structure  of  matter.  In  this  approach, 
antisymmetry  restrictions  are  enforced  subsequent  to  con¬ 
struction  of  the  many-electron  Hamiltonian  matrix  in  a  com¬ 
plete  orthonormal  spectral-product  basis.  Transformation  to 
a  permutation-symmetry  representation  obtained  from  the 
eigenstates  of  the  aggregate  electron  antisymmetrizer  is  seen 
to  enforce  the  requirements  of  the  Pauli  principle  ex  post 
facto,  and  to  eliminate  the  unphysical  (non-Pauli)  states 
spanned  by  the  product  representation.  Results  identical 
with  conventional  use  of  prior  antisymmetrization  of 
conhgurational  state  functions  are  obtained  in  applications 
to  many-electron  atoms.  The  development  provides  certain 
advantages  over  conventional  methods  for  polyatomic  mol¬ 
ecules,  and,  in  particular,  facilitates  incorporation  of  frag¬ 
ment  information  in  the  form  of  Hermitian  matrix  represen¬ 
tatives  of  atomic  and  diatomic  operators  which  include  the 
non-local  effects  of  overall  electron  antisymmetry.  An  exact 
atomic-pair  expression  is  obtained  in  this  way  for  polyatomic 
Hamiltonian  matrices  which  avoids  the  ambiguities  of  previ¬ 
ously  described  semi-empirical  fragment-based  methods  for 
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electronic  structure  calculations.  Illustrative  applications  to 
the  well-known  low-lying  doublet  states  of  the  H3  molecule 
in  a  minimal-basis-set  demonstrate  that  the  eigensurfaces  of 
the  antisymmetrizer  can  anticipate  the  structures  of  the  more 
familiar  energy  surfaces,  including  the  seams  of  intersec¬ 
tion  common  in  high-symmetry  molecular  geometries.  The 
calculated  H3  energy  surfaces  are  found  to  be  in  good  agree¬ 
ment  with  corresponding  valence-bond  results  which  include 
all  three-center  terms,  and  are  in  general  accord  with 
accurate  values  obtained  employing  conventional  high-level 
computational-chemistry  procedures.  By  avoiding  the 
repeated  evaluations  of  the  many-centered  one-  and  two- 
electron  integrals  required  in  construction  of  polyatomic 
Hamiltonian  matrices  in  the  antisymmetric  basis  states  com¬ 
monly  employed  in  conventional  calculations,  and  by  per¬ 
forming  the  required  atomic  and  atomic -pair  calculations 
once  and  for  all,  the  spectral-product  approach  may  provide 
an  alternative  potentially  efficient  ab  initio  formalism  suit¬ 
able  for  computational  studies  of  adiabatic  potential  energy 
surfaces  more  generally. 

Keywords  Potential  energy  surfaces  ■  Ab  initio 
calculations  •  Electronic  states  ■  Spectral  methods  ■ 
Antisymmetry  constraints 

1  Introduction 

Considerable  progress  has  been  reported  over  the  past 
four  decades  in  ab  initio  studies  of  the  complex  (Born- 
Oppenheimer)  potential  energy  surfaces  which  describe  the 
ground  and  low-lying  excited  electronic  states  of  molecules 
[1].  The  contributions  of  Professor  Mark  S.  Gordon  to 
this  enterprise  have  been  remarkably  comprehensive,  as 
typihed  by  his  central  role  in  the  development,  wide-spread 
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distribution,  and  continuing  collaborative  refinement  of  the 
well-known  GAMESS  suite  of  computer  codes  [2].  Accu¬ 
rate  ab  initio  methods  provided  in  this  and  other  compu¬ 
tational  resources  for  potential  energy  surfaces  commonly 
entail  repeated  calculations  of  large  numbers  of  one-  and 
two-electron  multicentered  integrals  over  explicitly  antisym¬ 
metrized  basis  states  in  construction  of  a  many-electron 
Hamiltonian  matrix  for  a  range  of  atomic  spatial  arrange¬ 
ments,  followed  by  determinations  of  the  energies  and  asso¬ 
ciated  eigenfunctions  of  selected  electronic  states  at  each 
atomic  arrangement,  and  construction  of  the  expectation  val¬ 
ues  corresponding  to  physical  properties  of  interest.  New 
methods  for  these  purposes  devised  in  a  spirit  of  contin¬ 
uing  collaborative  improvements  would  clearly  be  welcome, 
particularly  if  they  could  avoid  repeated  calculations  of  the 
individual  electronic  integrals  and  related  quantities  gen¬ 
erally  required  in  the  construction  of  potential  energy  sur¬ 
faces,  or  possibly  circumvent  entirely  the  determinations  of 
total  molecular  energies  (and  their  associated  differencing 
problems)  in  favor  of  calculations  of  atomic  and  interaction 
energies,  while  still  proving  applicable  to  both  ground  and 
electronically  excited  states  on  a  common  basis. 

In  a  contribution  to  the  systematic  improvement  of  stan¬ 
dard  quantum-chemical  methods,  an  alternative  perspective 
is  provided  in  the  present  report  on  enforcement  of  the  anti¬ 
symmetry  requirement  on  proper  atomic  and  molecular  wave 
functions  [3-11],  which  requirement  is  commonly  satisfied 
by  the  aforementioned  prior  constraint  on  the  representa¬ 
tional  basis  states  employed  in  variational  calculations  [1, 
2].  In  the  present  approach,  rather,  Pauli’s  exclusion  princi¬ 
ple  is  ignored  at  the  outset,  and  antisymmetry  is  enforced 
subsequent  to  construction  of  the  Hamiltonian  matrix  in  a 
spectral-product  representational  basis  familiar  from  the  the¬ 
ory  of  long-range  interactions  [12].  When  applied  to  the 
electronic  structures  of  atoms  as  an  introductory  pedagog¬ 
ical  example,  the  development  is  seen  to  be  equivalent  to 
the  familiar  Slater  approach  adopting  prior  antisymmetry 
of  many-electron  configurational  states  functions  [13],  pro¬ 
viding  some  degree  of  confidence  in  the  alternative  method 
more  generally.  In  applications  to  polyatomic  molecules,  the 
approach  provides  a  number  of  advantages  over  conventional 
methods,  and,  in  particular,  accommodates  the  incorporation 
of  ab  initio  atomic  and  diatomic  information  in  polyatomic 
Hamiltonian  matrices  in  a  rigorous  manner  which  avoids  cer¬ 
tain  difficulties  encountered  in  previously  described  semi- 
empirical  fragment-based  approaches  [14,15]. 

The  present  method  is  ultimately  based  on  conventional 
variational  calculations  in  L?  representations  of  many- 
electron  states,  guaranteeing  its  convergence  when  closure  is 
achieved,  in  which  limit  distinctions  between  use  of 
simple-product  and  explicitly  antisymmetrized  aggregate 
basis  states  are  known  to  become  formally  inconsequential. 
As  presented  here,  the  theory  is  applicable  to  polyatomic 


molecules  which  dissociate  into  neutral  species  on  their  low¬ 
est  lying  potential  energy  surfaces,  but  which  can  other¬ 
wise  involve  arbitrary  admixtures  of  covalent,  ionic,  van  der 
Waals,  and  metallic  interatomic  interactions.  Applications  of 
the  approach  to  ionic  systems,  or  to  neutral  systems  which 
can  give  rise  to  ion-pair  states  in  asymptotic  separation  lim¬ 
its,  require  only  minor,  largely  self-evident,  modifications  of 
the  present  description. 

The  atomic  spectral-product  basis  employed  in  the  devel¬ 
opment  is  known  to  span  the  totally  antisymmetric  irreduc¬ 
ible  representation  of  the  aggregate  electron  symmetric  group 
once  and  only  once,  making  it  suitable  for  analytical  and 
computational  studies  of  the  electronic  structures  of  many- 
electron  systems  [16-20].  A  matrix  representative  of  the  an- 
tisymmetrizer  constructed  in  the  spectral-product  basis  is 
employed  in  separating  the  totally  antisymmetric  and  non- 
totally-antisymmetric  (non-Pauli)  states  spanned  by  the  basis, 
and  in  correspondingly  isolating  the  physically  significant 
block  of  the  Hamiltonian  matrix  by  unitary  transformation. 
Hermitian  matrix  representations  of  atomic  and  atomic-pair 
operators  are  devised  in  this  way  for  polyatomic  molecules 
which  individually  have  well-defined  asymptotic  separation 
limits,  and  which  need  be  constructed  only  once  for  multi¬ 
ple  applications.  These  attributes  of  the  approach  facilitate 
development  of  an  ab  initio  unitary  transformation  formalism 
which  provides  the  aforementioned  exact  atomic -pair  rep¬ 
resentation  of  polyatomic  Hamiltonian  matrices,  including 
particularly  the  important  modihcations  of  the  bare  atomic- 
pair  interactions  due  to  the  aggregate  electron  antisymmetry 
required  upon  incorporation  into  the  polyatomic  system. 

The  general  theory  is  presented  in  Sect.  2,  where  the 
spectral-product  representation  of  electrons  is  dehned  and 
the  unitary  transformation  formalism  for  isolating  the  phys¬ 
ical  portion  of  the  Hamiltonian  is  described  for  both  atoms 
and  molecules.  Strategies  for  computational  implementation 
of  the  method  are  addressed  in  Sect.  3,  including  in  particular 
description  of  the  atomic-pair  implementation  and  illustra¬ 
tive  calculations  of  the  low-lying  doublet  electronic  states 
of  the  H3  molecule.  Concluding  remarks  summarizing  the 
salient  features  of  the  approach  are  given  in  Sect.  4. 

2  Theory 

The  spectral-product  approach  to  the  electronic  structures  of 
atoms  is  described  in  Sect.  2.1  as  a  pedagogical  introduc¬ 
tion  to  the  formalism,  whereas  its  application  to  polyatomic 
molecules  is  presented  in  Sect.  2.2. 

2.1  Spectral-product  formalism  for  atomic  structure 

Solutions  of  the  Schrodinger  equation  [21] 

//(r)^(r)  =  «'(r)-E  (1) 
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for  an  n -electron  atom  with  Hamiltonian  operator  [22] 


can  be  constructed  employing  a  formally  complete  square- 
integrable  (L^)  representation  in  the  outer-product  form  [16] 

$(r)  =  [</)(i)  (g)  (/)(2)  (g)  •  ■  ■  </>(n)}o  .  (3) 

The  row  vector  4>{i)  employed  in  Eq.  (3)  comprises  a 
denumerable  finite  or  infinite  basis  set  of  orthonormal  one- 
electron  spin-orbitals  specified  by  the  usual  quantum  num¬ 
bers  (n,  /,  mi,  s,  ms)  [22],  where  i  refers  to  the  spin  and  space 
coordinates  of  the  /th  electron  and  r  =  (i,  2,  . . . ,  n)  repre¬ 
sents  collectively  the  coordinates  of  all  n  atomic  electrons. 
The  particular  choice  of  spin-orbitals  is  irrelvant  for  the  pres¬ 
ent  development,  so  long  as  the  outer-product  (0)  Hilbert 
space  of  Eq.  (3)  can  include  a  suitable  closure  in  the  limit  of 
a  complete  spectral  basis  [23].  The  subscript  “O”  in  Eq.  (3) 
indicates  the  adoption  of  a  particular  ordering  convention  for 
the  sequence  of  product  functions  in  the  row  vector  $(r), 
the  consequences  of  which  convention  are  indicated  when 
appropriate  in  the  sequel. 

The  Hamiltonian  matrix  representative  of  the  operator  of 
Eq.  (2)  in  the  basis  of  Eq.  (3), 

H  =  ($(r)|H(r)|$(r)> 

=  X  (4) 

1=1  j=i+i 

includes  kinetic  =  (^>(r)|  —  (/j^/2m)V?|$(r)),  poten¬ 
tial  =  ($(r)|Ze^/r; |$(r)),  and  electron-interaction 
=  {$(r)|e^/r/y|$(r))  Hermitian  matrix  representa¬ 
tives  of  the  corresponding  operators.  The  dimensions  of  these 
matrices  are  determined  by  that  of  the  many-electron  basis 
of  Eq.  (3),  although  only  standard  one  and  two-electron  inte¬ 
grals,  and  their  products  with  the  unit  matrices  expressing 
the  orthonormality  of  the  remaining  n  —  1  and  n  —  2  orbitals, 
respectively,  are  required  in  evaluations  of  the  Hamiltonian 
matrix  of  Eq.  (4).  The  specific  forms  of  these  matrices  are 
determined  by  the  particular  ordering  convention  adopted  in 
Eq.  (3)  [16]. 

The  matrix  Schrodinger  equation  [2 1  ] 

H  Uh  =  Uh  E  (5) 

provides  approximate  or  exact  representations  of  the  row  vec- 
torofeigenfunctionsofEq.  (l)intheform^(r)  =  ^(r)UH, 
as  well  as  the  diagonal  matrix  of  associated  eigenvalues  E. 
The  energy  eigenspectrum  so  obtained  includes  both  physi¬ 
cally  significant  (totally  antisymmetric)  and  unphysical  (non- 
totally  antisymmetric  or  non-Pauli)  solutions,  all  of  which 
are  spanned  by  the  representation  of  Eq.  (3)  [17].  Accord¬ 
ingly,  this  particular  representation,  and  the  corresponding 


Hamiltonian  matrix  and  Schrodinger  equation  of  Eqs.  (4) 
and  (5),  can  encompass  descriptions  of  the  electrons  as  either 
Bose-Einstein,  Fermi-Dirac,  or  possibly  other  forms  of  iden¬ 
tical  particles,  necessitating  some  means  of  distinguishing 
among  these  possibilities  [24]. 

The  physical  and  non-Pauli  eigenstates  obtained  from 
Eqs.  (l)-(5)  can  be  separated  by  constructing  a  transformed 
Hamiltonian  matrix  which  has  a  physical  block,  a  non- 
Pauli  block,  and  vanishing  off-diagonal  blocks.  This  sepa¬ 
ration  is  accomplished  by  constructing  the  eigenstates  of  the 
n-electron  antisymmetrizer  in  the  representation  of  Eq.  (3) 
and  transforming  the  Hamiltonian  matrix  of  Eq.  (4)  to  this 
new  representation.  The  required  representation  matrix  is 

Pa  =  {0(r)|PA(«)l^(r’)),  (6) 

where  the  n-electron  Hermitian  antisymmetrizer  employed 
[24], 

n\ 

P^(n)^J^(-l)^»Pp,  (7) 

p=l 

is  left  unnormalized.  Alternatively,  adoption  of  wave- 
function  or  idempotent  normalization  in  Eq.  (7),  which  con¬ 
ventions  require  the  familiar  prefactors  ( 1  /  n ! )  ^  /^  and  ( 1  /  n ! ) , 
respectively,  does  not  materially  affect  the  ensuing  devel¬ 
opment,  except  for  multiplicative  changes  of  the  spectral 
interval  in  which  the  eigenvalues  of  the  antisymmetrizer  are 
placed. 

The  matrix  of  Eq.  (6)  is  found  to  comprise  a  series  of 
individual  non-zero  n\-hy-n  \  blocks  on  the  diagonal,  the  ele¬ 
ments  of  which  are  all  either  -t- 1  or  —  1 ,  and  a  remaining  set 
of  vanishing  blocks.  The  dimensions,  numbers,  and  natures 
of  the  non-zero  blocks  of  the  Pa  matrix  are  deteremined  by 
the  number  of  distinct  spin-orbital  products  employed  in  the 
representation  of  Eq.  (3),  and  by  the  presence  of  all  («!)  per¬ 
mutations  of  electrons  among  these  products.  This  structure 
ensures  that  any  finite  atomic  representation  in  the  form  of 
Eq.  (3)  provides  an  invariant  subspace  such  that  the  operation 
of  the  antisymmetrizer  is  closed  in  this  subspace 

pA{n)^(r)  =  $(r)  ■  Pa.  (8) 

As  a  consequence,  the  eigenfunctions  of  the  matrix  repre¬ 
sentative  of  the  antisymmetrizer  in  any  such  finite  space  will 
transform  correctly  as  either  totally  antisymmetric  or  non- 
totally  antisymmetric  basis  states  under  all  electron  permu¬ 
tations  [24]. 

The  eigenvalues  and  eigenfunctions  of  the  antisymmet¬ 
rizer  are  obtained  in  the  usual  way  from  the  unitary  transfor¬ 
mation  matrix  Up  that  diagonalizes  the  matrix  of  Eq.  (6), 

4  •  Pa  ■  Up  =  ,  (9) 

\  ^up  ^Ull / 
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where  the  physical  block  nllpp  contains  the  non-zero  ei¬ 
genvalues  and  the  unphysical  0„„  and  off-diagonal  blocks 
vanish  identically.  The  individual  non-zero  n!-by-n!  blocks 
of  Pa  indicated  above  are  each  diagonalized  separately  and 
contribute  a  single  non- vanishing  eigenvalue  (n!)  each  to  the 
physical  block  of  Eq.  (9).  All  other  eigenvalues  of  Pa  are 
identically  zero,  and  correspond  to  the  non-totally-antisym- 
metric  states  spanned  by  the  product  representation  of  Eq.  (3). 

The  eigenfunctions  of  the  antisymmetrizer  take  the  corre¬ 
spondingly  blocked  form 

«I>p(r)  ^  4»(r)  ■  Up  =  {{<^p(r)}p,  {<I»p(r)}„},  (10) 

consequent  of  the  ordering  of  the  columns  of  the  Up  matrix 
implied  by  Eq.  (9),  where  the  physical  states  satisfy 

PA(«){^p(r)}p  =  n!{‘I»p(r')}p  (11) 

and  the  non-Pauli  states  provide  the  null  results 
PA(n){^p(f )}/:>  =  0.  It  is  found  by  construction  that  the 
physical  eigenstates  of  Eqs.  (10)  and  (11)  are  linearly  inde¬ 
pendent  Slater  determinants  made  up  of  the  products  of  spin- 
orbitals  of  Eq.  (3)  which  can  provide  totally  antisymmetric 
states,  whereas  the  non-Pauli  eigenstates  are  totally  symmet¬ 
ric,  or  are  linear  combinations  of  other  (degenerate)  represen¬ 
tations  of  the  nth-degree  symmetric  group  [24].  The  number 
of  such  spin-orbital  products  leading  to  non-vanishing  Slater 
determinants  times  n!  gives  the  dimension  of  the  physical 
block  Ipp  appearing  in  Eq.  (9). 

Equations  (6)-(ll)  indicate  that  the  spectrum  of  PaW 
acting  in  the  domain  of  the  spectral-product  basis  of  Eq.  (3) 
is  that  of  a  compact  operator  [25],  the  value  zero  providing 
a  lower  limiting  point  of  accumulation  of  the  spectrum  of 
eigenvalues,  with  [$p(r  :  R)}u  the  associated  eigenstates. 
The  totally  antisymmetric  states  [$p(r  ;  R)}p  and  associ¬ 
ated  eigenvalues  n !  correspond  to  the  upper  limiting  point  of 
the  allowable  spectral  interval  (0,  n!). 

The  unitary  matrix  of  Eq.  (9)  provides  a  transformed 
Hamiltonian  matrix  in  the  form 

Hp^U;.H.Up=(*“P*^'"  ),  (12) 

V  »up  [HpImm/ 

where  [Hpjpp  and  [Hp}j,„  are  the  physical  (pp)  and  unphys¬ 
ical  (mm)  blocks  which  provide  the  totally  antisymmetric 
and  the  non-Pauli  solutions,  respectively,  and  the  off- 
diagonal  blocks  vanish  identically  in  accordance  with  Eqs.  (8) 
and  (10).  The  physical  Hamiltonian  matrix  of  Eq.  (12)  is 
identical  with  that  obtained  from  the  linearly  independent 
Slater  determinants  constructed  by  selecting  only  ordered 
spin-orbital  configurations  from  the  basis  of  all  such  prod¬ 
ucts  of  Eq.  (3).  Accordingly,  the  transformation  matrix  Up 
of  Eq.  (9)  obtained  from  the  matrix  representative  of  the 
antisymmetrizer  of  Eq.  (6)  incorporates  the  effects  of  elec¬ 
tron  permutation  symmetry,  in  the  absence  of  prior  basis-set 


antisymmetry  in  Eq.  (3),  in  the  Hamiltonian  matrix  of  Eq.  (4) 
through  the  transformation  of  Eq.  (12). 

2.2  Spectral-product  formalism  for  polyatomic  molecules 

The  foregoing  spectral-product  formalism  for  many-electron 
atoms  also  applies  with  some  important  differences  and  mod- 
ihcations  to  the  adiabatic  (Born-Oppenheimer)  electronic 
structures  of  polyatomic  molecules.  Specihcally,  in  the 
molecular  case  all  quantities  given  in  Sect.  2. 1  depend  explic¬ 
itly  upon  the  spatial  arrangement  of  atoms  in  the  molecule, 
specified  by  the  vector  R  =  (i?i,  R2, . . . ,  Rn),  where  Ra 
gives  the  position  of  the  atom  a.  In  order  to  follow  closely 
the  atomic  development  of  Sect.  2.1,  and  to  avoid  unnec¬ 
essary  notational  complexity,  the  molecular  development  is 
presented  here  for  a  collection  of  n  one-electron  (hydro¬ 
gen)  atoms  [17],  with  generalization  to  many-electron  atoms 
requiring  only  minor  additional  elaboration  [20]. 

The  spectral-product  representation  of  Eq.  (3)  for  poly¬ 
atomic  H„  molecules  becomes 

$(r  :  R)  =  |0^'^(i)  0  0®(2)  0  •  ■  ■  0^"^(n)|^  ,  (13) 

where  the  indvidual  spectral  states  0*^“^(i)  are  the  one- 
electron  spin-orbital  row  vectors  of  Eq.  (3),  with  the  elec¬ 
tron  coordinates  i  now  measured  relative  to  the  different 
atomic  origins  Ra .  As  in  the  case  of  Eq.  (3),  the  basis  of 
Eq.  (13)  spans  the  totally  antisymmetric  representation  of 
the  aggregate  n -electron  permutation  group  once  and  only 
once,  as  well  as  all  other  (non-Pauli)  representations  of  the 
nth-degree  symmetric  group,  in  the  limit  of  closure  [17]. 

The  molecular  Hamiltonian  operator  for  a  collection  of  n 
interacting  hydrogen  atoms  can  be  written  in  the  form  [cf., 
Eq.  (2)] 

n  n 

a=l  ^=01+1 


where  electron  i  (j)  is  arbitrarily  assigned  to  the  atomic  cen¬ 
ter  at  Ra  {Rp),  employing  common  notational  conventions 
[21].  A  decomposition  of  the  aggregate  Hamiltonian  opera¬ 
tor  similar  to  that  of  Eq.  (14)  is  also  possible  for  interacting 
many-electron  atoms,  although  sums  over  the  electrons  i  and 
j  assigned  to  the  atoms  a  and  P  are,  of  course,  required  in 
this  case  [16]. 
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The  matrix  representative  of  the  operator  of  Eq.  (14)  con¬ 
structed  in  the  basis  of  Eq.  (13)  [cf.,  Eq.  (4)], 

H(JS)  =  {$(r  :  R)\H(r  :  R)\^{r  :  R)) 

=  ^  V(“'^)(H„^)},  (15) 

is  seen  to  include  atomic  terms  =  (0(r  :  ii)|f/^“^(i)| 
0(r  ;  R))  and  atomic-pair  interaction  terms  V*^“’^'(i2a!yg)  = 
($(r  ;  R)\V^“’^\i\  j  :  Rap)\^(r  :  R))  which  provide 
Hermitian  matrix  representatives  of  the  corresponding  atomic 
and  interaction  operators.  The  atomic  energy  matrix 
is  independent  of  atomic  position  and  includes  only  diago¬ 
nal  one-electron  energies  when  the  orthonormal  atomic  spin- 
orbitals  are  chosen  to  be  atomic  pseudo-eigenstates, 

whereas  the  interaction  matrix  (Rafi)  is  generally  non¬ 
diagonal  and  depends  explicitly  upon  the  vector  separation 
Rap  of  the  two  indicated  atoms,  but  not  upon  the  individual 
laboratory-frame  positions  of  the  two  atoms,  nor  upon  the 
position  vectors  of  the  other  (n  —  2)  atoms  in  the  molecule. 
Evaluation  of  the  particularly  simple  Hamiltonian  matrix  of 
Eq.  (15)  evidently  requires  only  one-center  atomic  and  two- 
center  Coulombic  interaction  terms,  which  can  be  determined 
once  and  for  all  and  retained  for  repeated  applications.  These 
favorable  circumstances  are  consequences  of  the  orthonor¬ 
mality  of  the  spectral-product  basis  employed  in  the  absence 
of  prior  enforcement  of  overall  electron  antisymmetry,  the 
adoption  of  atomic  pseudo-states  in  the  representation,  and 
the  atomic  pairwise-additive  nature  of  the  interaction  terms 
in  the  Hamiltonian  operator  of  Eq.  (14). 

Following  the  development  of  Sect.  2.1,  the  Hermitian 
matrix  representative  of  the  antisymmetrizer  Pa  (n)  of  Eq.  (7) 
constructed  in  the  basis  of  Eq.  (13), 

Pa(H)  ^  {«^(r  ;  i2)|P,4(«)l‘l>(r  :  i2)> 

=  in\r^{PA(n)<S>ir  :  R)\PA(n)^(r  :  R)),  (16) 

is  seen  to  depend  explictly  on  the  atomic  arrangement  R, 
and  to  be  proportional  to  the  so-called  metric  matrix  of  the 
explicitly  antisymmetrized  product  basis  [21].  In  contrast  to 
the  atomic  matrix  of  Eq.  (6),  the  matrix  of  Eq.  (16)  is  a 
potentially  complicated  function  of  R,  requiring  detailed 
computational  evaluation.  Moreover,  the  operation  of  the  an¬ 
tisymmetrizer  is  not  closed  in  any  finite  representation  of  the 
form  of  Eq.  (13),  and  requires  a  complete  spectral  represen¬ 
tation  to  achive  closure  as  a  limiting  process  PA{n)^(r  :  R) 
^(r  :  R)  ■  PA(i2)  [cf.,  Eq.  (8)].  This  circumstance  is 
a  consequence  of  the  fact  that,  in  contrast  to  the  situation 
in  the  atomic  case  of  Eq.  (3),  the  representation  of  Eq.  (13) 
does  not  include  explicitly  all  n!  permutations  of  electron 
assignments  to  the  indicated  spin-orbital  product  states.  That 
is,  in  contrast  to  standard  treatments  [10,24],  in  which  “sym¬ 
metrical”  or  “democratic”  representations  of  all  electrons  are 


employed,  the  representation  of  Eq.  (13)  employs  different 
spin-orbital  basis  sets  to  describe  the  different  electrons  arbi¬ 
trarily  assigned  to  each  atom,  since  these  basis  sets  are  cen¬ 
tered  at  different  atomic  positions  Ra .  Accordingly,  electron 
permutations  are  included  in  Eq.  (13)  only  implictly  in  the 
limit  of  a  complete  spin-orbital  representation. 

The  unitary  transformation  matrix  Up  (R)  required  to  par¬ 
tition  the  Hamiltonian  matrix  of  Eq.  (15)  into  physical  and 
non-Pauli  blocks  is  obtained  from  the  diagonalization  [cf., 
Eq.  (9)] 

Up(fi)'  PA(fi)  Up(i2) 

_  t{Pd{R)}pp  OpM  (nWpp  OpjA 

“  V  V  {^a{R)]uu)  V  ■ 

Here,  the  eigenvalues  of  the  matrix  Pa(H)  are  partitioned 
into  an  upper  diagonal  block  [Pd(i2)}pp  containing  the  larg¬ 
est  eigenvalues  which  tend  to  the  upper  limiting  point  n\ 
of  the  spectrum  in  the  closure  limit,  and  a  lower  diagonal 
block  (Pd(i?)}uM  of  eigenvalues  which  tend  to  the  point  of 
accumulation  at  zero  in  this  limit  [17-20].  Similarly,  the 
permutation-symmetry-adapted  eigenstates  of  the  antisym¬ 
metrizer  are  obtained  in  this  case  as  a  limit  in  the  form  [cf., 
Eq.  (10)] 


$p(r  :  R)  =  $(r  ;  fi)  Up(i^) 

^  {(4»p(r:fi)}p,($p(r-:H)}„},  (18) 

where($p(r  ;  i2)}p  contains  the  totally  antisymmetric  states 
corresponding  to  the  non-zero  eigenvalues  [(Pd(J^)}/!p  — ^ 
nWpp]  of  the  antisymmetrizer,  and  (^p(r  ;  i?)}„  contains 
the  non-Pauli  states  corresponding  to  the  zero  eigenvalues 
[[Pd(i2)}Hii  ^  0„„]  of  the  antisymmetrizer,  in  accordance 
with  Eq.  (17).  Of  course,  the  partionings  of  Eqs.  (17)  and 
(18)  can  be  somewhat  subjective  in  any  finite  representation, 
in  contrast  to  the  situation  in  Eqs.  (9)-(12),  where  invariant 
physical  and  non-Pauli  subspaces  are  obtained  even  in  finite 
atomic  representations,  consequent  of  Eq.  (8). 

Finally,  the  unitary  matrix  Up(i2)  of  Eq.  (17)  is  employed 
in  constructing  the  transformed  Hamiltonian  matrix 

Hp(fi)  =  Up(i^)^  HiR)  Vp(R) 

a=l 

^  /(Hp(H)}pp  {lip(R)}pu\ 

\{lip(R)Up  [Hp(H)}„J 

/{Hp(i2)}pp  Opu  \ 

V  Oup  {lip(R)}uu)’  ^  ^ 

where  Hp“^(li)  and  Yp’^\R)  refer  to  transformations  of 
the  individual  atomic  and  interaction  matrices  in  Eq.  (15), 
(Hp(ii)}pp  and  (Hp(ii)}„„  are  the  physical  (pp)  and 
unphysical  (mm)  blocks  of  the  transformed  Hamiltonian 
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matrix  which  provide  the  antisymmetric  and  the  non-Pauli 
solutions,  respectively,  and  the  off-diagonal  blocks  vanish  in 
the  limit  (^)  of  spectral  closure  [17-20].  The  transformation 
of  Eq.  ( 1 9)  is  seen  to  incorporate  the  non-local  effects  of  over¬ 
all  electron  antisymmetry  in  the  atomic  ^ 
and  interaction  ^  matrices,  which 

consequently  now  depend  on  the  positions  of  all  atoms  in  the 
molecule.  The  Hermitian  atomic  and  pair-interaction  matri¬ 
ces  of  Eq.  (19)  individually  have  appropriate  dissociation 
limits,  with  the  entire  Hamiltonian  matrix  approaching  the 
sum  of  the  individual  atomic  terms  of  Eq.  (15)  in  the  com¬ 
plete  dissociation  limit  R  oo,  consequent  of  the  absence 
of  electron  exchange  in  this  limit  [PA(i2  ^  oo)  —y  I]  in  the 
denumerable  representation  of  Eq.  (13)  [20]. 


3  Computational  implementation  and  application 

Aspects  of  computational  implementation  and  numerical 
application  of  the  spectral-product  formalism  made  to  date 
are  summarized  in  Sect.  3.1.  The  atomic -pair-based  variant  of 
the  approach  is  reported  in  Sect.  3.2,  and  illustrative  numer¬ 
ical  applications  of  the  new  method  to  the  low-lying  doublet 
states  of  the  H3  molecule  are  given  in  Sect.  3.3. 

3.1  Implementation  strategies 

The  Hamiltonian  matrices  of  Eqs.  (15)  and  (19)  related  by 
the  unitary  transformation  matrix  of  Eq.  (17)  provide  iden¬ 
tical  energy  eigenvalues  and  eigenfunctions  which  converge 
to  those  of  the  operator  of  Eq.  (14)  in  the  limit  of  closure 
of  the  spectral-product  representation  of  Eq.  (13).  A  num¬ 
ber  of  strategies  devised  for  verifying  the  ideas  and  theorems 
on  which  the  spectral-product  development  is  based,  and  for 
determining  the  desired  Hamiltonian  eigenspectra,  has  been 
employed  to  date  in  selected  cases  [16-20].  Brief  descrip¬ 
tions  of  progress  made  in  these  approaches  aid  in  under¬ 
standing  the  issues  which  must  be  overcome  in  devising  a 
generally  applicable  computational  implementation  of  the 
spectral-product  formalism. 

Although  computations  based  directly  on  Eq.  (15),  which 
entail  only  evaluations  of  the  indicated  atomic  and  atomic- 
pair  interaction  terms,  appear  attractive,  previously  described 
technical  complications  associated  with  the  presence  of 
unphysical  continua,  spuriously  degenerate  solutions,  and 
slow  convergence  in  long-range  separation  limits  mitigate 
against  adopting  such  an  approach  in  general  [20].  In  the 
special  case  of  two-electron  systems,  however,  the  electron 
spin  functions  can  be  factored  out,  the  remaining  allowable 
spatial  wave  functions  are  either  symmetric  or  antisymmetric 
under  transposition  of  the  electron  coordinates,  and  unphys¬ 
ical  (non-Pauli)  or  spuriously  degenerate  solutions  are  not 


present.  Accordingly,  computational  studies  of  the  H2  mole¬ 
cule  based  on  Eqs.  (13)-(15)  can  be  performed  in  appropri¬ 
ately  chosen  orbital-product  basis  sets,  the  results  of  which 
exhibit  monotone  convergence  with  basis  set  expansion  to 
the  well-known  lowest  lying  singlet  and  triplet  energies  [18]. 
Examination  of  the  electronic  charge  densities  of  these  states 
clarifies  the  manner  in  which  the  spectral-product  represen¬ 
tation  can  describe  the  accumulation  or  depletion  of  negative 
charge  between  the  atoms  relative  to  the  undisturbed  atomic 
charge  distributions,  providing  thereby  chemical  bonding  or 
anti-bonding  in  the  two  states  entirely  in  the  absence  of  com¬ 
monly  employed  classical  chemical  structures  to  describe  the 
wave  functions  [18].  Additionally,  the  expectation  values  of 
the  electron  transposition  operator  P,  /  employing  the  singlet 
and  triplet  spectral-product  wave  functions  are  seen  to  con¬ 
verge  smoothly  with  basis  set  expansion  to  the  appropriate 
values  (±  1)  in  each  case  [18].  These  results  provide  confi¬ 
dence  in  the  development  of  Eqs.  (13)-(15)  in  applications 
to  the  prototypically  important  electron  pair  bond. 

Evaluations  of  the  Pa(P)  matrix  of  Eq.  (16)  for  the  H2 
molecule  in  appropriately  chosen  basis  sets  over  the  range 
R  <  10  ao  of  atomic  separation  verify  the  convergence  of  the 
eigenvalues  of  Eq.  (17)  to  the  upper  (n\  ^  2)  and  lower  (0) 
end  points  of  the  allowed  spectral  interval  in  this  case  [19, 
20] .  Construction  of  the  transformed  Hamiltonian  matrix  of 
Eq.  (19)  verifies  the  blocking  indicated  there,  corresponding 
to  separation  of  the  singlet  and  triplet  states  in  the  absence  of 
unphysical  representations  in  this  two-electron  case.  Eurther- 
more,  the  rate  at  which  the  off-diagonal  blocks  of  the  Ham¬ 
iltonian  matrix  which  potentially  connect  the  two  different 
spin  symmetries  become  negligible  with  increasing  basis  is 
made  quantitative  in  the  calculations  [19,20]. 

As  an  alternative  to  constructing  the  unitary  transforma¬ 
tion  matrix  of  Eq.  ( 17)  and  the  Hamiltonian  matrix  of  Eq.  (19) 
for  the  H2  molecule,  the  representation  of  Eq.  (13)  is  aug¬ 
mented  with  explicitly  antisymmetric  Heitler-London  singlet 
or  triplet  functions  and  the  Hamiltonian  matrix  of  Eq.  (15) 
correspondingly  augmented  with  an  additional  row  and  col¬ 
umn  [19,20].  This  implementation  is  equivalent  to  the  devel¬ 
opment  of  Eqs.  (16)-(19),  but  entails  performing  separate 
calculations  for  the  singlet  and  triplet  states  in  a  manner 
which  serves  to  accelerate  the  convergence  of  each  individual 
calculation.  In  this  way,  accurate  potential  energy  curves  are 
obtained  for  the  H2  molecule  which  converge  smoothly  from 
their  Heitler-London  approximations  to  the  essentially  exact 
Kolos-Wolniewicz  curves  [26],  providing  further  confidence 
in  the  development  of  Eqs.  (16)-(19)  [19,20]. 

Applications  to  weakly  bound  (van  der  Waals)  atomic 
clusters  provide  useful  first  steps  in  implementation  of  the 
spectral-product  method  for  many-electron  atoms  more  gen¬ 
erally.  In  the  case  of  inert-gas  clusters  doped  with  a 
single  optically  active  atomic  radical,  attention  is  focused 
on  efficient  construction  of  the  individual  Coulombic 
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atomic-radical/inert-gas  pair-interaction  matrices  required  in 
Eq.  (15)  [16].  In  this  approach,  similar  in  spirit  to  atoms- 
and  diatomics-in-molecules  and  related  methods  [14, 15]  but 
differing  significantly  in  its  implementation  [16],  a  unitary 
transformation  from  diatomic  to  spectral-product  represen¬ 
tations  is  employed  to  decompose  the  adiabatic  diatomic 
states  and  potential  energy  curves  calculated  employing  stan¬ 
dard  methods  [1,2]  into  diabatic-like  atomic-product  states 
and  the  associated  interaction  energy  matrices  of  Eq.  (15). 
Although  successfully  employed  in  “on-the-fly”  construc¬ 
tions  of  ground  and  electronically  excited  potential  energy 
surfaces  in  Monte  Carlo  and  molecular  dynamics  simula¬ 
tions  of  singly  doped  inert-gas  clusters  [27],  the  absence  of 
the  transformation  of  Eqs.  (16)  to  (19)  in  this  approach  for 
van  der  Waals  clusters  renders  it  unsuitable  for  more  gen¬ 
eral  applications  involving,  in  particular,  covalent  and  ionic 
bonding. 

3.2  Atomic  pair  implementation 

The  particular  form  of  the  Hamiltonian  matrix  of  Eq.  (19) 
suggests  the  possibility  of  systematic  evaluation  of  the  indi¬ 
vidual  atomic  and  interaction-energy  terms  appearing  in 
Eq.  (1 5),  for  which  evaluations  can  be  performed  once  and  for 
all  independent  of  any  particular  atomic  aggregate  arrange¬ 
ment,  followed  by  construction  of  the  aggregate-specific 
transformation  matrix  Up  (R)  of  Eq.  ( 1 7).  It  is  convenient  for 
both  pedagogical  and  notational  simplicity  to  describe  this 
approach  for  three  interacting  hydrogen  atoms  following  the 
notation  of  Sect.  2.2,  with  the  more  general  many-electron- 
atom  case  reported  separately  elsewhere  [28]. 

In  the  three-atom  case,  the  basis  of  Eq.  (13)  takes  the  form 

^(r  :R)  =  {0(“Hi)®0(''^(2)®0(")(3)j^,  (20) 

where  R  =  (Ra,  Rb,  Rc)  and  r  =  (i,  2,  3)  are  the  atomic 
and  electronic  coordinates,  respectively.  The  corresponding 
Hamiltonian  matrix  of  Eq.  (19)  includes  three  atomic  terms 
and  three  atomic-pair  interaction  terms  of  the  forms 

H^,“^(i^)  =  Up(i^)^■H(“^■Up(i^)  (21) 

vj,“'^^(i2)  =  Up(i2)^  ■  ■  Up(H),  (22) 

where  the  atomic  Hamiltonians  in  Eq.  (21)  (a  =  a,  b, 
or  c)  are  evaluated  employing  standard  atomic  spectral  cal¬ 
culations  [22]. 

In  order  to  incorporate  specifically  diatomic  calculations 
in  evaluation  of  the  atomic-pair  interaction  energy  matrices 
Vp“’^^(i2)  of  Eq.  (22)  (afi  =  ab,  ac,  or  be),  the  transfor¬ 
mation  matrix  Up(i?)  of  Eq.  (17)  appearing  there  is  fac¬ 
tored  [=  U^’^^(i2ayg)  ■  (i?)]  into  two  parts,  where  the 

unitary  matrix  Up“’^\i?a;^)  diagonalizes  the  matrix  repre¬ 
sentative  P^^’^\Rap)  of  the  two-electron  antisymmetrizer 


constructed  in  the  spectral-product  basis  of  Eq.  (20), 

PA’^\Raii)  ^  mr  :  R)\P^/-^\2)mr  :  R)) 

=  (23) 

and  the  unitary  matrix  Up^“^\i?)  described  further  below 
completes  the  transformation  from  spectral-product  to  three- 
electron  states  of  good  permutation  symmetry  [Eq.  (18)]. 
Evaluation  of  the  matrix  P^^'^\Rap)  defined  in  Eq.  (23) 
involves  forming  the  outer  product  of  the  two-electron  met¬ 
ric  matrix  (Rap)  with  the  unit  matrix  arising  from 
the  orthogonality  of  the  third  “by-stander”  spin-orbital  row 
vector,  the  elements  of  the  resulting  full  matrix  {Rap) 

arranged  in  accordance  with  the  ordering  convention  implied 
by  the  subscript  “O”  in  Eq.  (20)  [16]. 

The  transformation  matrix  \Jp’^\Rap)  that  diagonalizes 
P^^’^\Rap)  of  Eq.  (23)  provides  an  intermediate  “diatomic” 
basis  set  in  the  form  [cf.,  Eq.  (18)] 

«I>^"’^^(r  ;  R) 

^  «^(r  ;  R)  ■  U^“'^^(i2„/i) 

=  {{{0<“)(i)  ®  0(^>(j)}  ■  u^^’^\Rap)}  ®  <l>^^'>{k)}^ 

^  {(<’^^(r  :  R)}p^^,  [<’^^(r  :  i^)}„„^},  (24) 

where  Up’^\Rap)  refers  to  the  strictly  diatomic  portion  of 
the  transformation  matrix  Up‘^'^^(iiQ.^)  =  (Up“’^^(i?Q.^)  (g) 
the  physical  states  (^p“’^\r  :  R)]p^p  are  antisym¬ 
metric,  the  unphysical  states  :  R)}u^p  are  symmet¬ 

ric  in  the  electron  coordinates  i  and  j  in  the  closure  limit, 
and  the  subscript  pap  (uap)  designates  the  dimension  of  the 
physical  (unphysical)  diatomic  subspace.  Note  that  designa¬ 
tion  of  the  basis  states  of  Eq.  (24)  as  “diatomic”  refers  to  the 
specifically  diatomic  calculations  required  in  their  construc¬ 
tion  employing  only  the  atomic-pair  row  vector  {<j)^‘^\i)  (g 
which  diatomic  states  are  subsequently  augmented 
by  multiplication  with  the  third  by-stander  spin-orbital  row 
vector  as  indicated  explicitly  in  the  second  line  of 

Eq.  (24)  [16]. 

Employing  the  foregoing  factoring  Up  (i?)  =  Up“’^^(i2Q;^) 
■Up^“^^  (R),  the  pairwise-atomic  interaction-energy  matrix  of 
Eq.  (22)  can  now  be  written  in  the  form 

V^,“’^^(i^)  =  lJp‘‘''\R)^  ■  \^p'^\Rap)  ■  Vp‘‘^\R),  (25) 

where  the  transformed  interaction-energy  matrix 

^^p'^\Rap)  ^  V^f^\Rap)^  ■  V(“’«(i2„^)  ■  V^p^^\Rap) 

=  {<'^^(r  :  i2)|y(“'^)(i;  j  :  i^„^)|<•^^(r  :  i^)) 

(26) 
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is  evaluated  in  the  indicated  permutation-symmetry-adapted 
diatomic  basis  of  Eq.  (24),  and  should  not  be  confused  with 
the  matrix  Vp“’^^(i?)  of  Eq.  (22)  and  the  left-hand  side  of 
Eq.  (25).  Noting  that  the  operator  j  :  Rap)  can 

be  written  as  the  difference  of  two  diatomic  Hamiltonian 
operators  ji  :  Rap)  —  :  Rap  -»  oo)] 

which  are  symmetric  functions  of  the  electron  coordinates  i 
and  j,  the  matrix  of  Eq.  (26)  is  seen  to  be  block  diagonal  in 
the  diatomic  physical  (pap)  and  unphysical  (uap)  labels,  in 
accordance  with  Eq.  (24). 

The  matrix  required  to  complete  the  transfor¬ 

mation  of  Eqs.  (25)  and  (26)  is  obtained  from  diagonalization 
of  the  matrix  representative 

^  {<“’^>(r-  :  H)|P,4(3)|<“’^^(r  :  R)) 

=  V^p’^\Rap)^  ■  Pa(-R)  ■  V^p'^\Rap)  (27) 

of  the  three-electron  antisymmetrizer  PaO)  constructed  in 
the  diatomic  basis  of  Eq.  (24),  where  Pa(H)  is  its  represen¬ 
tation  in  the  spectral-product  basis  of  Eq.  (20).  Accordingly, 
it  is  clear  from  this  expression  that  if  Up^“^\il)  diagonalizes 
P^^°‘^\R)  of  Eq.  (27),  the  product  matrix  lJ^^\Rap) 

■  Up^“^\i2)  diagonalizes  PA(i2). 

The  three  sets  of  three-electron  permutation-symmetry- 
adapted  states  obtained  from  the  foregoing  development  [cf., 
Eq.  (18)], 

:  R)  =  :  R)  ■  Vp‘‘^\R) 

=  4>(r  :  R)  ■  V^p^\Rap)  ■  V^i^^^kR) 

^  :  R)}p,  :  fi)}4, 

(28) 

differ  only  in  the  sequence  of  permutation  symmetry 
adaption  employed  in  their  construction.  Accordingly,  the 
functions  of  Eq.  (28)  are  all  equal  to  each  other  and  to  the 
defining  states  of  Eq.  (18)  ;  R)  $p(r  ;  ii)] 

in  the  limit  of  closure  in  Eqs.  (24)  and  (27).  In  this  limit,  the 
three  matrices  Up“^\i?)  can  be  conveniently  obtained  in 
the  form  =  Up“^^(ii„^)^  ■  Up(il),  where  the  dia¬ 

tomic  matrices  Up“^\i?Q,^)  for  aP  =  ab,  ac,  or  be  can  each 
be  constructed  from  an  antisymmetrizer  matrix  P^^^\Rap) 
which  depends  only  upon  the  scalar  separation  Rap  employ¬ 
ing  standard  angular-momentum  considerations  to  provide 
the  dependences  upon  the  orientation  of  the  vector  separa¬ 
tions  Rap  [16]. 

Einally,  since  the  antisymmetrizer  matrix  of  Eq.  (27)  is 
block  diagonal  in  the  diatomic  labels  pap  and  Uap  in  the 
closure  limit,  the  off-diagonal  blocks  {V^°'^\R)} ua^ pap 

of  the  matrix  vanish  identically 


in  this  limit.  Accordingly,  the  physical  block  of  the 
interaction-energy  matrix  of  Eq.  (25)  takes  the  form 

{V^.“’^’(H)}pp  =  ■  {v<.“’^^(fi„^)) 

^  PPa^  ^  h’a^Pa^ 

■  {Up^“'’\fi)l  ,  (29) 

I  i  Pa^P 

where  {'Vp'^hRap)} pappap  ^^e  physical  block  of  the 
interaction-enegy  matrix  of  Eq.  (26),  and  the  the  transfor¬ 
mation  matrix  V^p°‘^\R)} p^p p  is  generally  rectangular  in 
view  of  the  generally  different  dimensions  of  the  physical 
subspaces  pap  and  p  (<  pap)  of  Eqs.  (24)  and  (28),  respec¬ 
tively.  Equation  (29)  indicates  that  unphysical  diatomic  states 
do  not  contribute  to  the  interaction  terms,  and  that  only  the 
physical  block  {Up^“^^  (R))pc,pp  transformation  matrix 

Up^“^^(i?)  is  required  in  the  development. 

A  final  expression  for  the  physical  block  of  the 
Hamiltonian  matrix  of  Eq.  (19)  is  obtained  from  the  fore¬ 
going  exact-atomic-pair  development  of  Eqs.  (20)-(29)  in 
the  form 

n  n 

{Hp(R)}pp^J^  {{<^(i2)}pp-f  XfVp’^^(i2)}pp), 

(30) 

where  the  atomic  terms  are  evaluated  using  Eq.  (21)  and  the 
pair-interaction  terms  are  given  by  Eq.  (29)  and  the  associ¬ 
ated  preceeding  expressions. 

3.3  Illustrative  computational  applications 

Computational  applications  of  the  atomic -pair  formalism 
reported  here  are  based  on  the  previously  demonstrated 
general  equivalence  between  the  physical  subspace  of  a  spec¬ 
tral-product  representation  and  the  linearly  independent 
canonically  orthogonalized  subspace  of  the  associated  explic¬ 
itly  antisymmetrized  form  of  such  a  representation  [20].  In  a 
simple  model  illustration  adopting  this  approach,  the  lowest 
lying  doublet  potential  energy  surfaces  of  the  H3  molecule  are 
studied  here  in  a  minimal-basis-set  representation  [29,30]. 

The  physical  subspace  of  the  diatomic  states  of  Eq.  (24), 
although  formally  defined  in  terms  of  the  spectral-product 
representation  of  Eq.  (20)  employed  there,  is  equivalently 
obtained  in  the  present  computational  application  from 
explicitly  antisymmetric  basis  states  in  the  familiar  canon¬ 
ically  orthogonalized  form  [20,21] 

{4»^“'^^(r  :  i2)l  ^ 

^  J  Pa^ 

^{P^«./5)(2)  o(r  :  R)  -  U^.“■^^(i^„^)}p„, 
■l^^:'^\R.p)rpllip,  (31) 
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where  the  transformation  and  eigenvalue  matrices  employed 
are  obtained  from  the  antisymmetrizer  matrix 
of  Eq.  (23)  given  in  the  preceeding  development.  In  turn, 
Eq.  (31)  is  employed  in  evaluation  of  the  physical 
blocks  and  of  the 

interaction-energy  and  antisymmetrizer  matrices  of  Eqs.  (26) 
and  (27),  respectively,  both  of  which  are  required  in  con¬ 
struction  of  the  final  expressions  of  Eqs.  (29)  and  (30)  for  the 
atomic-pair  form  of  the  Hamiltonian  matrix  . 

An  expression  similar  to  that  of  Eq.  (31)  can,  of  course, 
also  be  employed  for  the  physical  subspace  of  the  aggregate 
states  of  Eq.  (18)  for  H3  in  the  form 

{$p(r  ;  R)}p 

^{.Pa(3)  $(r  :  R)  ■  \]v{R)]p  ■  ,  (32) 


The  unitary  transformation  matrix  of  Eq.  (16)  in  the  rep¬ 
resentation  of  Eq.  (33)  is 

Pa(H) 

!  1  ^ab^ac^bc  ^ab^ac^bc 

—  ^ab^ac^bc  1  ^ac  ^ab^ac^bc  ^ 

^ab^ac^bc  ^ac  ^ab^ac^bc  1  ^ 

(34) 

where  =  s{Ra^')  =  (ljQ;(i)|H^(i))  <  1  is  the  one- 
electron  li  orbital  overlap  integral.  The  eigenvalue  Pd(J?) 
and  transformation  Up(i?)  matrices  of  Eq.  (17)  are  obtained 
from  diagonalization  of  the  matrix  of  Eq.  (34)  following  the 
development  of  Eqs.  (16)-(18)  in  the  forms 


where  the  indicated  transformation  and  eigenvalue  matrices 
are  those  of  Eq.  (17).  Employing  this  expression  in  evaluat¬ 
ing  the  Hamiltonian  matrix  in  the  usual  variational  approach 
gives  rise  to  the  multicenter  integrals  [29]  the  atomic-pair 
development  is  designed  to  avoid.  That  is,  in  contrast  to 
Eq.  (32),  the  atomic-pair  expression  for  antisymmetric 
aggregate  basis  states  is  given  by  Eq.  (28),  in  which  only 
the  diatomic  states  indicated  there  are  made  explicitly  anti¬ 
symmetric  in  accordance  with  Eq.  (31),  with  the  remaining 
antisymmetry  in  these  functions  accounted  for  in  the  limit 
by  the  indicated  transformation  matrix 
Accordingly,  it  is  seen  that  the  atomic-pair  implementation 
of  the  spectral-product  approach  can  be  expected  to  converge 
in  the  limit  to  results  obtained  from  use  of  the  valence-bond- 
like  expression  of  Eq.  (32),  also  taken  to  a  suitable  limit  of 
convergence  [30]. 

Adopting  the  foregoing  approach  in  a  minimal-basis-set 
description  of  the  H3  molecule,  and  restricting  attention  to 
Ms  =  -1-1/2  states,  the  spectral-product  basis  employed  in 
Eq.  (31)  takes  the  three-term  form 


«I>(r  :  R)  =  l^,(l)Hft(2)H,(3)[a(l)a(2)/l(3), 


{Pd(.R)}l,2  —  1  ^ab^bc^ca  i 


+(s^.b-4cf  +  (sl-4cf} 


2  x2i1/2 


{Pd(i2)}3  —  ^  -  Sab 


Sbc  4”  ^SabSbcSac^ 


(35a) 

(35b) 


and 


Us(i^) 

/ ^  cos  0(R)  -  ^  sin  0(R)  ^  sin  0(i?)  +  ^  cos  0(R) 

_  ^cos&(R)  +  ^sin&(R)  ^  sin0(i2)  —  ^  cos0(ii)  ^ 

^  -^cos0(fi)  sin  0(H) 

(36a) 


where 


0(i?)  =  -  tan 


-1 


(sL  -  He  +  sic) 


(36b) 


a(l)/l(2)a(3),/l(l)a(2)a(3)},  (33) 

where  a  (i)  andj6(i)  are  the  Pauli  spin  functions,  and  = 

1  ia  ( I  ■f’i — R-a  I )  refers  to  a  H  orbital  centered  on  atom  a .  In  this 
approximation,  only  the  three  indicated  spin  couplings  con¬ 
tribute  to  the  expressions  of  Eq.  (20)-(30),  with  the  diatomic 
interactions  treated  at  the  Heitler-London  level  of  approx¬ 
imation  spanning  only  the  lowest  lying  singlet  and  triplet 
states. 


The  two  eigenvalues  given  by  Eq.  (35a)  and  the  first  two 
columns  of  the  Up(H)  matrix  correspond  to  doublet  states 
of  H3,  whereas  Eq.  (35b)  and  the  last  column  of  Up(H)  cor¬ 
respond  to  a  quartet  state,  the  functions  of  Eqs.  (28)  and  (32) 
proving  5'=l/2,  3/2;  Ms  —  +1/2  spin  multiplets  in  this  sim¬ 
ple  representation. 

The  corresponding  diatomic  antisymetrizer 
and  transformation  Up“’^^(7?a^)  matrices  of  Eqs.  (23)  and 
(24),  respectively,  are  obtained  from  appropriate  limits  of 
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Eqs.  (34)-(36a),  (36b)  in  the  forms 


where  the  individual  atomic-pair  interaction  terms  have  been 
combined  into  a  total  interaction-energy  matrix  of  the  form 


„{a,b) 


{Rab)  = 


\it^\Rab)  = 


^t''\Rac)  = 


0  0 
1  -^Ib 


I  0 

VI 


\b 

1  0\ 


0 

0 

1 

-5^ 

■'ac 

-^2 

■'ac 

0 

0 

0 

1 

0 

0 

1 

^/2/ 

1 

<^2 

^bc 

-4c 

1 

0 

0  1 

1 

0 

J_\ 

72 

V2\ 

1 

0 

1 

V2 

V2 

0 

1 

0/ 

(37a) 


(37b) 


(37c) 


and  the  eigenvalues  of  the  three  diatomic  P^^^\Rap)  matri¬ 
ces  all  take  the  identical  form 


0  0  \ 

0  l-.2^  0 

0  0  1-4/ 


(38) 


Construction  of  the  diatomic  states  of  Eq.  (3 1 )  employing  the 
Vp’^\Rap)  matrices  of  Eqs.  (37a)-(37c)andthePj“’^\i?„^) 
matrix  of  Eq.  (38)  shows  them  to  be  products  of  antisymmet¬ 
ric  Heitler-London  diatomic  states  {ab,  ac,  or  be)  of  *  E+ 
or  symmetry  with  a  by-stander  li'  atomic  spin  orbital 
(c,  b,  or  a)  of  appropriate  electron  spin. 

In  accordance  with  the  discussion  following  Eq.  (28),  the 
transformation  matrices  required  to  evaluate  Eq.  (29)  are 
obtained  from  the  matrices  of  Eqs.  (36a),  (36b)  and  (37a)- 
(37c),  which  information  completes  evaluation  of  the 
Hamiltonian  matrix  of  Eq.  (30)  for  the  H3  molecule  in  a 
minimal-basis  representation.  In  this  case,  the  Hamiltonian 
can  be  written  in  a  form  that  emphasizes  the  separate  con¬ 
tributions  from  atomic  energies  and  atomic-pair  interaction 
potentials. 


Hs(i2)  =  3Eul  +  VsiR)^  ■  \(R)  ■  VsiR),  (39) 


/3 


\(R)^ 


Vab  +  vie"’  -H  V, 


,(av) 


V, 


(diff) 


be 

Miff) 


yldiff) 

fv!f+ 

yidiff) 

^ab 


,(av) 


^Vbc 


yidiff) 

ac 

yidiff) 

''ab 

hvr+ 


/(av) 


(40) 


Here,  ^Vap  =  ^V{Rap)  and  ^Vq,^  =  ^ViRap)  are  the 
triplet  and  singlet  Heitler-London  interaction-energy  curves, 
respectively,  and  =  [^V(Rap)  +^V{Rap)]/2  is  the 

average  and  =  [^y(I?«yg)—^y(/?o!;3)]/2  one-half  the 

difference  of  these,  all  of  which  potentials  vanish  in  the  limit 
Rap  ^ 

It  is  seen  that  the  effect  of  the  transformation  matrix  Us  (R) 
in  Eq.  (39)  is  to  modify  the  individual  pairwise  potentials 
of  Eq.  (40)  to  include  non-pairwise-additive  or  three-body 
interactions  in  the  Hamiltonian  matrix,  in  common  with  ear¬ 
lier  semi-empirical  attempts  to  incorporate  three-center  terms 
in  pair-interaction-based  approximations  to  the  H3  molecule 
energy  surface  [3 1-36].  Detailed  comparisons  with  these  ear¬ 
lier  approaches  are  of  considerable  interest,  and  are  reported 
separately  elsewhere  [28].  Carrying  out  the  indicated  matrix 
multiplications  in  Eq.  (39)  shows  that  {Hs(li)}i3  = 
{Hs(i?)}23  =  0,  as  expected  of  Hamiltonian  matrix  elements 
connecting  doublet  and  quartet  states.  Additionally,  the  ele¬ 
ment  {Hs(i2)}i2  is  generally  small,  and  vanishes  identically 
in  high-symmetry  atomic  arrangements  [Cjv,  Dooh,  D3/,]. 
Of  course,  when  the  entire  matrix  of  Eq.  (39)  is  retained,  its 
eigenvalues  are  identical  with  those  obtained  in  the  absence 
of  the  transformation  matrix  Us(i2),  whereas  more  generally 
only  a  small  physical  block  of  the  transformed  Hamiltonian 
is  employed  [Eq.  (19)]. 

In  the  spirit  of  the  blocking  of  Eq.  (19),  quantitative  com¬ 
parisons  are  provided  here  of  the  predictions  of  the  two 
doublet-state  diagonal  elements  of  the  model  Hamiltonian 
matrix  of  Eqs.  (39)  and  (40)  with  valence-bond  calcula¬ 
tions  which  include  explicit  three-center  terms  obtained  from 
the  representation  of  Eq.  (32)  [30].  The  calculations  per¬ 
formed  employ  a  12-term  Gaussian  representation  of  the  H 
atomic  hydrogen  orbital  which  provides  an  orbital  energy 
(—0.4999  au)  in  good  agreement  with  the  correct  value,  as 
well  as  explicitly  antisymmetric  singlet  and  triplet  molecular 
hydrogen  eigenstates  and  a  ground-state  equilibrium  inter¬ 
atomic  separation  {Rg  =  1.642ao)  and  total  energy 

(—1.1157  au)  in  good  agreement  with  the  accepted  Heitler- 
London  values  [21]. 

Figure  1  shows  contour  plots  of  the  two  doublet  eigen¬ 
values  of  the  Pa  (ii)  matrix  of  Eq.  (34)  obtained  in  the  li^  H3 
representation  for  non-symmetric  collinear  atomic  arrange¬ 
ments.  Although  the  eigenvalues  of  the  Pa(7?)  matrix  will 
converge  in  accordance  with  Eq.  (17)  in  arbitrarily  large  basis 
sets,  the  minimal-basis  results  of  Fig.  1  are  nevertheless  seen 
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Fig.  1  Constant-value  contours  for  the  two  doublet  eigenvalues  of  the 
Pa(.R)  matrix  ofEq.  (34)  for  collinear  atomic  configurations  of  H3  con¬ 
structed  in  a  representation,  employing  an  increment/decrement  of 
0.05  between  adjacent  contours.  The  coordinate  axes  give  the  distances 
Rab  and  Rac  between  the  central  atom  located  at  the  origin  and  the 
outer  atoms,  measured  in  atomic  units  (uq)-  In  a  the  contours  increase 
uniformly  in  value  from  that  labeled  1.05,  except  as  R^b  =  Rac  0 
inside  the  saddle  point,  in  which  region  the  values  decrease  monotoni- 
cally  to  zero,  whereas  in  b  the  contour  values  decrease  uniformly  from 
that  labeled  0.95 

to  provide  qualitatively  useful  information.  Specifically,  the 
eigenvalue  surface  of  Fig.  la  has  a  maximum  along  the  line 
R^l,  —  at  1.80  UQ,  decreases  rapidly  to  zero  in  the  limit 
R^l,  —  R^^  — >•  0,  approaches  1  in  the  three-body  break-up 
limit  Rah  —  Rac  00,  and  goes  to  2  in  the  diatomic  limits 
Rab  0,  Rac  — ^  00  and  Rab  — >•  00,  Rac  0.  The  smaller 
doublet  eigenvalue  surface  shown  in  Fig.  1(b)  is  evidently 
less  structured  than  that  of  Fig.  1  (a)  and  vanishes  in  both  the 
diatomic  and  triatomic  united-atom  limits.  The  vanishing  of 
the  eigenvalues  of  Fig.  1  as  Rab  and  Rac  -»  0  herald  the  for¬ 
mation  of  non-Pauli  states  in  the  representation  of  Eq.  (33)  or 
of  linearly  dependent  doublet  states  in  the  representation  of 
Eq.  (32),  indicating  in  either  event  that  the  minimal  1a^  basis 
does  not  support  physically  significant  states  in  the  united- 
atom  limit  H3  ^  Li.  The  Pa  (i?) -matrix  eigenvalue  surfaces 
evidently  identify  atomic  configurations  in  which  improper 
states  arise,  provide  information  complementary  to  the  asso- 


Fig.  2  Constant-energy  contours  for  the  ground-  and  first-excited 
doublet-state  potential  energy  surfaces  of  H3  for  non-symmetric  collin¬ 
ear  atomic  arrangements,  obtained  from  Eqs.  (39)  and  (40)  constructed 
in  a  li^  representation,  employing  an  energy  decrement/increment  of 
0.5  eV  between  contours.  The  coordinate  axes  are  as  in  Fig.  1,  and 
the  zero  of  energy  employed  (— 1.50au)  is  that  of  the  three-atom  dis¬ 
sociation  limit  (H  -r  H  -r  H).  a  depicts  the  familiar  saddle  and  dia¬ 
tomic  structures  of  the  ground-state  energy  surface  in  H3  [33],  whereas 
the  excited-state  energy  surface  of  b  is  montonically  increasing  with 
decreasing  separations  Rab  and/or  Rac 

elated  energy  surfaces,  and  can  anticipate  the  presence  of 
crossings  in  the  corresponding  potential  energy  surfaces,  as 
demonstrated  in  further  detail  below. 

Eigure  2  provides  constant-energy  contour  plots  of  the 
two  lowest  lying  doublet  energy  surfaces  in  H3  obtained 
from  Eqs.  (39)  and  (40)  for  the  non-symmetric  collinear 
atomic  arrangements  of  Eig.  1 .  These  well-known  energy  sur¬ 
faces  are  presented  here  to  illustrate  the  extent  of  agreement 
obtained  with  corresponding  conventional  valence-bond  cal¬ 
culations  in  a  basis  [29,30],  and  to  demonstrate  their 
general  accordance  with  accurate  quantum-chemistry  cal¬ 
culations  [37-39].  Specifically,  the  calculated  saddle-point 
energy  (—1.574  au)  of  Eig.  2a  compares  well  with  that 
obtained  from  valence-bond  calculations  (— 1.583au)  in  the 
same  basis,  and  the  H3  saddle  point  geometry  {Rab  —  Rac  — 
2.18flo)  is  in  general  accord  with  the  accepted  position 
{Rab  —  Rac  —  1.757ao)  [37-39].  The  present  calculation 
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places  the  saddle  point  approximately  0.0735  an  (2.00 eV) 
below  the  three-body  (H  -i-  H  -i-  H)  breakup  energy  of  —  1 .50au, 
and  thus  provides  approximately  50%  of  the  accepted  accu¬ 
rate  saddle-point  binding  energy  of  0.1570  au  (4.273  eV) 
relative  to  this  limit  [37-39].  The  calculated  saddle-point 
energy  relative  to  the  two-body  break-up  energy  (H  -f  H2) 
provides  a  1.14eV  barrier  to  the  exchange  reaction  which 
is  in  accord  with  the  valence-bond  calculations  (0.876  eV) 
but  is  about  three  times  larger  than  the  accepted  value  of 
0.4166eV  [37-39],  emphasizing  the  well-known  sensitivity 
of  this  energy  difference  to  the  computational  approximation 
employed  [33]. 

The  structure  of  the  ground-state  energy  surface  of  Fig.  2a 
is  evidently  complementary  to  that  of  the  associated  PA(i2)- 
matrix  eigenvalue  surface  of  Fig.  la,  with  larger  values  of 
the  PA(i?)-matrix  eigenvalues  corresponding  to  generally 
lower  energy  values.  Similarly,  the  energy  surface  for  the 
excited  doublet  state  shown  in  Fig.  2b  is  seen  to  be  monoton- 
ically  increasing  from  the  three-atom  break-up  limit  as  the 
diatomic  and  triatomic  united-atom  limits  are  approached, 
complementary  to  the  monotonically  decreasing  form  of  the 
associated  PA(i2)-matrix  eigenvalue  surface  of  Fig.  lb.  The 
failure  of  the  li^  representation  to  provide  an  accurate 
ground-state  H3  energy  surface  near  the  saddle  point  can  be 
associated  with  the  absence  of  physically  significant  states 
near  the  united-atom  limit,  as  anticipated  by  the  antisymmet- 
rizer  eigenvalue  surfaces  of  Fig.  1 . 

In  Fig.  3  are  shown  quantitative  comparisons  of  the  first 
two  doublet  energies  in  FI3  with  corresponding  three-center 
valence-bond  calculations  for  symmetric  collinear  (Rab  = 
Rac)  atomic  arrangements.  Evidently,  the  results  of  Eqs.  (39) 
and  (40)  for  the  ground-state  energy  surface  are  in  good 
agreement  with  the  Is^  valence-bond  calculations  which 
include  only  covalent  structures  (o),  and  are  also  in  accord 
with  valence-bond  results  which  include  all  possible  ionic 
configurations  (•)  in  the  basis  [30].  The  first  excited  dou¬ 
blet  state  in  H3  is  seen  to  be  repulsive,  in  good  agreement 
with  the  valence-bond  calculations  and  in  qualitative  accord 
with  previously  reported  accurate  calculations,  although  this 
excited  doublet  is  known  to  undergo  an  avoided  crossing 
with  a  higher  lying  state  which  is  not  included  in  the  pres¬ 
ent  development  [40].  Although  the  results  of  Eig.  3  for  the 
ground-state  energy  surface  appear  to  bound  the  valence- 
bond  results  from  above,  this  is  only  a  happenstance.  Indeed, 
similar  apparent  boundedness  is  not  found  for  the  excited 
doublet  state  in  Eig.  3,  although  both  results  are  in  satisfac¬ 
tory  agreement  with  the  valence-bond  calculations. 

Figure  4  depicts  contour  plots  of  the  two  doublet-state 
Pa  (i?) -matrix  eigenvalues  in  H3  for  C2v  arrangements  in 
which  the  Rat  =  Rac  separations  are  co-varied  and  the  apex 
angle  9  is  varied  from  30°  to  180°.  A  seam  of  intersection 
of  these  two  surfaces  indicated  by  the  cusps  in  the  contours 
is  evident  in  the  figures  at  9  —  60°,  which  corresponds  to 


Fig.  3  Ground  and  first-excited  doublet  state  energies  in  H3  for  sym¬ 
metric  collinear  atomic  arrangements.  The  solid  curves  depict  the  values 
obtained  from  Eqs.  (39)  and  (40)  constructed  in  a  representation, 
whereas  the  data  points  are  obtained  from  conventional  valence-bond 
calculations  in  a  basis  set  including  {filled  circles)  and  excluding 
(open  circles)  ionic  configurations  [30].  The  zero  of  energy  employed 
(— l.SOau)  is  that  of  the  three-atom  dissociation  limit  (H  +  H  +  H),  as 
in  Fig.  2 


the  familiar  energy  degeneracy  associated  with  £>3/,  symme¬ 
try  [33].  The  dashed  and  solid  portions  of  the  contours  in 
the  two  figures  depict  the  continuous  diabatic-like  surfaces 
which  pass  smoothly  through  one  another  at  the  9  —  60°  line 
of  intersection,  in  accord  with  the  degeneracy  predicted  by 
Eq.  (35a).  In  panel  (a)  the  higher  value  eigensurface  (solid 
lines)  is  to  the  right  of  the  60°  line  and  the  lower  value  eigen¬ 
surface  (dashed  lines)  to  the  left  of  this  line  of  intersection. 
Correspondingly,  panel  (b)  shows  the  reversed  arrangement, 
with  the  lower  eigenvalue  surface  (dashed  lines)  to  the  right 
and  the  higher  (solid  lines)  to  the  left  of  the  60°  line. 

The  associated  energy  surfaces  obtained  from  Eqs.  (39) 
and  (40)  for  the  ground-  and  first-excited  doublet  states  of 
H3  for  C2v  atomic  arrangements,  shown  in  Fig.  5,  similarly 
depict  the  expected  D^h  degeneracy  at  9  —  60°  and  the 
associated  surface  crossing  along  the  cusps  in  the  contours, 
complementary  to  the  corresponding  seam  of  intersection 
anticipated  by  the  PA(i?)-matrix  eigenvalues  of  Fig.  4.  This 
surface  crossing  arises  in  the  present  representation  from  the 
aforementioned  vanishing  of  the  off-diagonal  Hamiltonian 
matrix  element  {Hs(i?)]i2  of  Eqs.  (39)  and  (40)  at  high  sym¬ 
metries. 


^  Springer 


Theor  Chem  Acc 


30  60  90  120  150  180 

0  (degrees) 


Fig.  4  Constant-value  contours  for  the  two  doublet-state  eigenvalues 
of  the  Pa{R)  matrix  of  Eq.  (34)  constructed  in  a  li^  representation 
for  C2v  atomic  arrangements  (Rab  =  Rac),  where  9  is  the  apex  angle 
of  the  triatomic  H3  configuration,  employing  an  increment/decrement 
of  0.1  between  adjacent  contours.  The  cusps  in  the  contours  of  both 
panels  indicate  the  presence  of  a  seam  of  intersection  at  6  =  60°  in  the 
two  continuous  surfaces  identified  by  the  solid  and  dashed  lines  corre¬ 
sponding  to  the  two  eigenvalues  of  Eq.  (35a),  with  the  contour  values 
decreasing  uniformly  from  right  to  left  in  a  and  decreasing  uniformly 
from  top  to  bottom  in  b,  as  discussed  further  in  the  text 


In  Fig.  6  are  shown  comparisons  of  the  two  lowest  lying 
doublet  energies  in  H3  with  coresponding  valence-bond  cal¬ 
culations  for  the  C2v  atomic  arrangements  of  Fig.  5  in  which 
Rab  and  Rac  are  held  fixed  at  the  ground-state  Heitler-London 
separation  of  1.64  qq  and  6  is  varied  from  20°  to  180°.  Evi¬ 
dently,  the  present  results  are  in  good  agreement  with  the  cor¬ 
responding  three-center  1j ^  valence-bond  calculations  which 
include  only  covalent  structures  (o).  By  contrast,  although  the 
present  results  are  also  in  accord  with  valence-bond  results 
which  include  all  possible  ionic  configurations  (•)  in  the  li^ 
basis  [30],  the  undulation  between  120°  and  180°  evident 
upon  inclusion  of  such  terms  in  the  valence-bond  calcula¬ 
tions  is  evidently  not  present  in  the  l^^  implementation  of 
the  atomic-pair  formalism.  The  crossing  of  the  two  surfaces 
when  9  —  60°  is  evident  in  the  figure,  and  in  accord  with  the 
valence-bond  calculations.  As  in  Fig.  3,  the  apparent  upper 
bound  provided  by  the  present  results  on  the  valence-bond 


0  (degrees) 

Fig.  5  Constant-energy  contours  for  the  two  lowest  lying  doublet  states 
of  H3  for  C2ti  atomic  arrangements  (Rah  =  Rac),  obtained  from 
Eqs.  (39)  and  (40)  constructed  in  a  li^  representation,  where  6  is 
the  apex  angle  of  the  triatomic  configuration,  employing  an  energy 
increment/decrement  of  0.5  eV  between  adjacent  contours.  The  zero 
of  energy  (—  1 .50au)  is  that  of  the  three-atom  dissociation  limit  (H  -r  H 
+  H).  The  cusps  in  the  contours  of  both  panels  indicate  the  familiar  £>3/, 
symmetry  seam  of  intersection  at  0  =  60°  in  the  two  continuous  energy 
surfaces  identified  by  the  solid  and  dashed  lines  [33],  with  the  contour 
energy  values  increasing  uniformly  from  top  to  bottom  in  both  panels, 
a  depicts  the  lowest  lying  surface  only  for  values  of  6  >  60° ,  whereas 
b  depicts  the  lowest  lying  surface  only  for  6  <  60°,  as  is  discussed 
further  in  the  text 

calculations  in  the  case  of  the  ground  state  in  Fig.  6  is  a  hap¬ 
penstance,  with  a  similar  apparent  upper  boundedness  not 
present  in  the  case  of  the  excited  doublet  state. 


4  Concluding  remarks 

The  theoretical  development  reported  here  provides  an  alter¬ 
native  perspective  on  electronic  structure  calculations  for 
atoms  and  molecules  and  other  forms  of  matter.  Antisym¬ 
metry  restrictions  are  enforced  in  this  approach  subsequent 
to  construction  of  the  many-electron  Hamiltonian  matrix  of 
an  atom  or  molecule  in  an  orthogonal  spectral-product  basis. 
Transformation  to  a  permutation-symmetry  representation 
obtained  from  the  eigenstates  of  the  aggregate  electron 
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Fig.  6  Ground  and  first-excited  doublet-state  energies  in  H3  for  C2v 
atomic  arrangements  (Rab  =  Rac  =  l-64ao),  where  6  is  the  apex 
angle  of  the  triatomic  configuration.  The  solid  curves  depict  the  values 
obtained  from  Eqs.  (39)  and  (40)  constructed  in  a  Ij^  representation, 
whereas  the  data  points  are  obtained  from  conventional  valence-bond 
calculations  in  a  Is^  basis  set  including  (filled  circles)  and  excluding 
(open  circles)  ionic  configurations  [30],  The  zero  of  energy  (— 1.50au) 
is  that  of  the  three-atom  dissociation  limit  (H  -I-  H  -r  H),  as  in  Fig.  5. 
The  crossing  point  of  the  two  curves  corresponds  to  D3;,  symmetry,  as 
discussed  further  in  the  text 


antisymmetrizer  enforces  the  requirements  of  the  Pauli  prin¬ 
ciple  in  this  approach.  Results  are  obtained  in  applications  to 
many-electron  atoms  which  are  identical  with  the  use  of  an¬ 
tisymmetrized  configurational  state  functions  in  variational 
calculations,  providing  some  degree  of  confidence  in  the 
soundness  of  the  method  more  generally.  In  applications  to 
the  polyatomic  molecules,  the  development  accommodates 
the  incorporation  of  fragment  information  in  the  form  of 
Hermitian  matrix  representatives  of  atomic  and  diatomic 
operators  which  include  the  non-local  effects  of  overall  elec¬ 
tron  antisymmetry  on  the  individual  atomic  and  pairwise- 
atomic  interaction  terms. 

Standard  computational  methods  employing  explicitly 
antisymmetric  diatomic  wave  functrions  are  employed  in 
constructing  the  pair-interaction  matrices  and  other  quan¬ 
tities  required  in  an  implementation  of  the  method  in  the 
pro  to  typically  important  H3  molecule.  The  eigensurfaces  of 
the  matrix  representative  of  the  electron  antisymmetrizer  are 
seen  to  identify  the  presence  of  non-Pauli  states  in  the 
representation  employed,  to  complement  the  structures  of 
the  more  familiar  energy  surfaces,  and  to  anticipate  the  pres¬ 


ence  of  seams  of  surface  intersection  associated  with  high- 
symmetry  molecular  geometries.  The  minimal-basis-set 
results  reported  provide  energy  surfaces  in  good  agreement 
with  corresponding  three-center  valence-bond  calculations, 
and  are  in  general  accord  with  the  results  of  accurate  quan¬ 
tum-chemistry  calculations. 

The  formalism  reported  avoids  repeated  calculations  of 
the  one-  and  two-electron  many-centered  integrals  generally 
required  in  construction  of  polyatomic  Hamiltonian  matrices 
employing  the  antisymmetric  basis  states  commonly  adopted 
in  calculations  of  potential  energy  surfaces.  Rather,  the 
approach  entails  atomic  and  molecular  calculations  which 
can  be  performed  once  and  for  all  and  retained  for  repeated 
applications,  in  combination  with  construction  of  an 
aggregate-specific  transformation  matrix  designed  to 
incorporate  the  non-local  interaction  effects  associated 
with  enforcement  of  complete  electron  antisymmetry. 
Accordingly,  the  present  development  can  possibly  provide 
an  alternative  ab  initio  approach  suitable  for  computational 
applications  to  polyatomic  molecules  more  generally. 
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